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Abstract 

We present a high statistics, quenched lattice calculation of the S-parameters Bb^ 
and Bb^, computed at lowest order in the HQET. The results were obtained using a 
sample of 600 quenched gauge field configurations, generated by Monte Carlo simula- 
tion at /3 = 6.0 on a 24^ x 40 lattice. For the light quarks the SW-Clover action was 
used; the propagator of the lattice HQET was also tree-level improved. Our best esti- 
mate of the renormalization scale independent S-parameter is I3bj = 1-03 it 0.06ib0.18. 
Bb^ has been obtained by using "boosted" perturbation theory to calculate the renor- 
malization constants which relate the matrix elements of the lattice operators to the 
corresponding amplitudes in the continuum. Due to the large statistics, the errors 
in the extraction of the matrix elements of the relevant bare operators are rather 
small. The main systematic error, corresponding to ±0.18 in the above result, comes 
from the uncertainty in the evaluation of the renormalization constants, for which the 
one-loop corrections are rather large. The non-perturbative evaluation of these con- 
stants will help to reduce the final error. We also obtain Bb^/Bb^ = 1-01 ib 0.01 and 
pBBBjflBB, = l.38±0m. 



1 Introduction 



The study of B^-B^ and B^-B^ mixings may narrow the constraints on the elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, which in turn serve to test the 
Standard Model, and to detect possible new physics beyond it [|| . The results of the 
phenomenological analyses which involve B^-B^ mixings |^ strongly depend on two 
non-perturbative inputs: the value of the decay constant of the i?-meson, /b, and 
the value of the renormalization group invariant i?-parameter, Bb- These parameters 
enter in the theoretical predictions for and x^, to be defined below, through the 

mixing parameter = fs \J Bb, £,b corresponding, up to trivial factors, to the matrix 
elements of the relevant AB = 2 operators. The i?-meson decay constant has not 
been measured yet. However, several theoretical evaluations of Jb using numerical 
simulations of lattice QCD and QCD sum rules |^ have been published to date. 
The S-parameter, which is equal to one in the vacuum saturation approximation, has 
also been estimated using the same non-perturbative methods ||2|, |3| . 

In this paper a high statistics lattice calculation of the i?-parameter, at lowest 
order in the HQET, is presented. The relevant matrix elements have been computed 
on a statistical sample of 600 gauge field configurations. This has allowed us to study 
two- and three-point correlation functions at large time-distances with relatively small 
statistical errors, thus reducing the systematic errors due to the contamination of 
higher excited states. For the bare four-fermion operators, it was possible to isolate a 
good signal by using cube and double-cube smeared axial currents as the interpolating 
fields for the 5-mesons, and to extract their matrix elements with small statistical 
errors. The major source of systematic error is presently given by the uncertainties 
in the evaluation of the renormalization constants of the different operators, which 
are known only at first order in perturbation theory. Indeed, as explained below, by 
renormalizing the operators in two different ways, which are equivalent up to O(a^), 
we obtain either 

Bb^ = 1.21 ± 0.08 Method 1 (1) 

or 

Bb^ = 0.86 ± 0.06 Method 2 . (2) 
By combining eqs. ([l|) and (|2|), we derive our final result for this quantity 

= 1.03 ±0.06 ±0.18, (3) 

where the last error is an estimate of the uncertainty due to the contribution of higher- 
order terms to the renormalization constants. 

The ratio of i?-parameters Bb^/Bb^, which may be obtained by studying the de- 
pendence of the i?-parameter on the light quark mass, has also been considered. In 

^ The subcripts d or s referring to the light quark flavour will be written explicitely only when necessary 
for the discussion. 
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this ratio, most of the uncertainties due to higher-order terms cancel and we quote 



1.01 ±0.01. (4) 



By using fss/fBa = 1.17ib0.03, as derived from the study of the two-point correlation 
functions on the same set of configurations, and Mb^/Mb^ = 1.017 we then obtain 



r,. = -f--f-^ = 1.43 ±0.07. (5) 



flMlBB, 
flMlBB, 

Tsd can also be obtained by taking directly the ratio of the matrix elements of the 
A.B = 2 operator, computed at the appropriate value of the quark masses. In this case 
we get 

fl Ml Bb, 

f^^M^^R =1.35 ±0.05. (6) 

Jb^^b.Bb, 

smaller than the result r^^ = 1.54(13) obtained using propagating Wilson quarks in- 
stead of the HQET [§ 0. The experimental values of Mb^ and Mb, were taken from 
ref. @. 

Only a non-perturbative calculation of the renormalization constants, which relate 
the relevant lattice operators to the continuum one, can reduce the present uncertainty 
in the prediction of the 5-parameters from the HQET. This computation is in progress 
and the results will be published elsewhere. 

The plan of the paper is the following. In sec. |2| the relevant formulae for the 
definition of Bb are given; in sec. |^ the procedure followed to compute the i?-parameter 
is presented; in sec. § the numerical calculations of the relevant matrix elements are 
described and the main results of this study are discussed; in the conclusions an outlook 
of possible developments of the calculations of this paper is given. 



2 Master formulae for B^—B^ mixings 

mixings are usually expressed as 

Xd,s = , (7) 

where (AM)^^_^ is the mass difference of the B'^^-B^ ,, system and ^b^s — ^/'^Bds 
the average total width. Using the effective /S.B = 2 Hamiltonian discussed in ref. 
one finds 

Xd,s = TB^^, ^ r]B Mb^^, fBa,s Soixt) \Vt{d,s)? , (8) 



^ In ref. they also presented a result which includes the error due to the extrapolation to a — > 0. 
Since our result has been obtained only at ~ 2 GeV, we compare it with the result of ref. [Q obtained 
by fitting Tsd to a constant in a. 
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where ijb ~ 0.55, So{xt) is a function of the ratio xt = mf/M^r, rrit and M\y are the 
top quark and W masses respectively, Gp the Fermi constant, Vt[d,s) tlie Cabibbo- 



Kobayashi-Maskawa matrix element and Csds = fs^s v^^ds- 

The renormalization group invariant parameter Bb^^ is defined by 



where 



t H 

47r 



^BdJ/^) , (9) 



fi is the operator renormalization scale and /^^ ^ is the -B^^s-meson decay constant 

(0|6(0)7^75g(0)|i?(p-)) = ip^fB . (11) 

Oi is the renormalized four-quark operator 

Ol = ( b{x) -f^ (1 - 75) q{x) ) ( 7^ (1 - 75) (?(rE) ) , (12) 

and b{x) and q{x) are the heavy and light quark fields respectively. 

The factor J5 can be written in terms of the one- and two-loop anomalous dimen- 
sions, 7^^^ and 7^^\ of the operator Ol] if OiifJ-) is renormalized in the MS'-scheme 
one gets 

7(0) = 4, ^{i) = _7+ln/, /3o = ll-^n/, 



38 7W /3i 

A = 102--n,, J„, = 2^-2_. ,13) 



/3o,i are the one- and two-loop coefficients of the f3 function and nj the number of 
active light quark flavours. Uf = 5 for /i > m^. Notice that, strictly speaking, Bb is 
renormalization scale independent only up to the next-to- leading order. 

Since is larger than the current values of the lattice cut-off, in order to predict 
physical quantities relevant in i?-physics a possibility is to use the lattice version of the 
HQET 1^, which allows a systematic expansion in inverse powers of the heavy quark 
mass. It this case, to obtain the physical matrix elements of the effective Hamiltonian, 
one has to match appropriate lattice bare operators to the operator Ol{^), which is 
renormalized in the continuum "full" theory. The matching procedure is conventionally 
splitted in two separate steps 0: 

i) the matching of the continuum HQET to the full theory; 

ii) the matching of the lattice HQET to its continuum counterpart. 



As explained in the first of refs. Q this is, however, unecessary. 
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The two-loop Next-to-Leading Order (NLO) anomalous dimension necessary for 
step i) has been computed using the MS dimensional regularization in ref. |^]. As 
noticed in ref. |^], the calculation was, however, not complete because some elements of 
the one-loop mixing matrix were missing. This calculation has been recently completed 
in refs. pO| , 11 1. The final result can be written as 

dL{m,) = C7i(/x2)(5L(/i) + Cs{^i^)^s{^J^) , (u) 



where OL{mh) is the operator in the full theory renormalized in the MS'-scheme at the 
a scale equal to 6-quark mass (OL("^b) = OiifJ- = rub)) and 



Ol 
Os 



( 7^ (1 - 75) q{x) ) ( h^~^ (x) 7^ (1-75) q{x) 
( h{x) (1 - 75) q{x) ) ( /i(-) {x) (1 - 75) q{x) 



(15) 



are the operators of the effective theory, renormalized at the scale fi. The fields h and 
h^^'^ create a heavy-quark or annihilate a heavy anti-quark state respectively. The 
coefficient functions are given by ||lO|, 0] 



asiml) 



d2 



.(0) 
721 



- (0) - (0) 

722 - Tn 



Cs{ml) 



d2 



Csim, 



where, in naive dimensional regularization 



b) > 



ClK) 



with 



1 - 14- 



di 



40) 



Csiml) 



J = /3] 



7ii 



2/3o "'/3o 2/3o 

The elements of the one-loop anomalous mixing matrix are given by 



7ii 



-(0) 
7l2 







"(0) 
721 



-(0) 
722 



and 



7!!' 



805 



52 . 64 

— Tr^-\ 

27 9 



(16) 



(17) 
(18) 

(19) 
(20) 
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We now consider step ii), i.e. the matching of the lattice to the continuum HQET. 
Owing to the breaking of chiral symmetry induced by the Wilson term, the discretized 
version of Ol, mixes with two new lattice operators as follows p^-[p!^] 



OL{^^) = + In(aV') + DL]j OUa) 

+ ^DnOnia) + ^DNON{a) (21) 
47r 47r 

Osifi) = (^l + ^[^ln(aV) + ^5]) 05(a) + ... (22) 
where Oi{a) denote a bare lattice operator; 

+ ( h{x) 7^ (1 + 75) q{x) ) ( h^-^ (x) 7^ (1 - 75) q{x) [ 
+ 2 ( hix) (1-75) q{x) ) ( h^-^ (x) (1 + 75) q{x) ) 
+ 2 {h{x){l+-f5)q{x)) [U-Hx){l--f5)q{x)) ; (23) 

is the lattice coupling constant and Dl,... ,Dr are constants which have been 
computed in one-loop lattice perturbation theory. The leading logarithmic correction 
to O5 in eq. (p^ ) is taken into account by renormalizing the operators at a scale = 1/a; 
the term Dg and the . . . represent finite terms of 0{as)- Since Os enters only at the 
NLO, they do not need to be computed, cf. eqs. (|l^) and ([TtI). Two remarks are in 
order here: 

• The value of the constants Di depend on the light quark action used in the 



numerical simulation. In our case the tree- level improved SW-Clover action |15] 



was used and the operators were improved by rotating the light-quark propagators 
In the calculation of the coefficients Dj, the heavy-quark field renormalization 



constant computed in ref. |17]-[19| was used. This is the definition which is 



consistent with improvement |Ty] and with the method used to extract the matrix 



elements of the operators Oj(a), see also the discussion in refs. [20, 

The details of the numerical calculation of the coefficients Di can be found in ref. |14]. 
Here, we only quote the final results: 

Dl = -22.5 Dr = -5.4 Dn = -14.0 (24) 
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Lattice perturbation theory behaves rather badly due to the presence of tadpole-like 
diagrams which are typical of the lattice case 0. For logarithmically divergent operators, 
two remedies to this problem have been proposed in the recent years, namely boosted 
perturbation theory l23] or the non-perturbative renormalization of the operators on 



quark and gluon states, in a fixed gauge [24|. Only the first method can be used. 



since the non-perturbative calculation of the renormalization constants relevant in the 
present study has not been performed yet. 

In eq. (|2lD, for comparison with the recent calculation of the UKQCD collaboration 
af' = 6/47r/3Mo has been used, uq is a measure of the average link variable, 
uq = {8Kc)^^, with Kc the value of the Wilson hopping parameter K at which the pion 
mass vanishes. In the SW-Clover case, the numerical value of the "boosted" coupling 
given above is very close to the boosted coupling defined in terms of the elementary 



plaquette (see ref. |23| for the different definitions). We have also used = av{q*), 
where av{q*) was also introduced in ref. [^, and q* is an appropriate scale, which 
can be extracted in one- loop perturbation theory. For q*, the value computed for the 



heavy-light axial-vector current for the Wilson action in ref. [25| was used, q*a = 2.18, 
with a-^ = 2 GeV, since the appropriate result for the present case is not known. 

Putting the corrections of steps i) and ii) together, see eqs. (|l^), (^), (p2|) and 
(p^), the renormalized operator OL{m\)) can be written in terms of the bare lattice 
operators Oj(a) as follows: 

OL(mfe) = Zo^ OL{a) + Zoj, Onia) + Zq^ On {a) + Zqs Os{a) . (25) 

The renormalization constants Zoi are obtained from products of the coefficients in 
eqs. ([TBI), (^) and (p^), which are computed independently and contain terms of 0{as)- 
Thus there are two different ways in which one can organize the final results, i.e. by 
including (Mi) or excluding (M2) next-to-next-to- leading terms of O(a^) (without log- 
arithms). Since the corrections are large, different choices will result in i?-parameters 
which differ by about 28%. With the present state of art this is an intrinsic systematic 
error that can be reduced only by incresing the order of the perturbative calculation 
in the continuum and applying the non-perturbative renormalization to the lattice 
operators. 

In eq. ([l^), the Zq^ have been evaluated by using vrib = 5 GeV, ~ ^^O 

MeV, /X = 2 GeV ~ 1/a and Uf = A. At the NLO this choice of the parameters gives 
as{'m\) = 0.1842, which is also the value used in eq. (P), since the matching between 
the full theory and the effective one is made at ^ = rrib- 



In eq. (21) three possible values for were used, and the mixing coefficients Zq^ 
were computed with both options Mi and M2 ■ The values of the Zq^ , together with 
the corresponding values of are given in table |l[ In the same table we also give the 
renormalization constant of the heavy- light axial current, computed in ref. |^6| at the 
NLO, which is needed for the calculation of the i?-parameter, see sees. ^ and ^ 
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The inaccuracy of lattice perturbation theory has been demonstrated explicitely in several cases where 



a non-perturbative determination of the renormalization constants was possible [021 G3| 
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Options 


Zol 


Zor 


Zon 


Zos 




Standard 
.0796 


Ml 


.7949 


-.0317 


-.0822 


-.1229 


.9035 


Ma 


.7651 


-.0394 


-.1021 


-.1229 


.8989 


ag Boosted uq 
.1458 


Ml 


.6850 


-.0581 


-.1506 


-.1229 


.7932 


Ma 


.6286 


-.0722 


-.1871 


-.1229 


.7847 


Boosted V 
.1800 


Ml 


.6283 


-.0717 


-.1859 


-.1229 


.7363 


Ma 


.5581 


-.0891 


-.2310 


-.1229 


.7257 



Table 1: Renormalization constants for different choices of the lattice coupling constants and 
options, see the text. 

Two observations are necessary at this point. 

• In ref. [^l[ they have used the option Ma, without, however, including the next- 
to-leading corrections of step i). Moreover, they took for Og the value given by 
its leading logarithmic expression. This is not consistent at the next-to-leading 
order and introduces a further (and easily avoidable) systematic effect of ~ 10% 
in the final result. 

• Since, in absence of a better control of perturbation theory, the uncertainties are 
so large, it is not worth at this point to vary the value of ^qcd^ or to worry 
about "quenching", or other more subtle questions. 



3 Computation of Bb on the lattice. 

In order to obtain Bb, we compute the following two- and three-point correlation 
functions 

C^R'it) =Y^{Q\A^{x,t)A^'\^,m , (26) 

X 

CBfit^M) = E (0K(xi,ti)O,(0,0)<'^(x2,ta)|0) , (27) 

where i =L,R,N,S denotes one of the operators on the r.h.s. of eq. (|25|). The labels 
R,R' =L,S,D correspond to different interpolation operators used for the i?-mesons 

^li^) = E^(^*)>^5^(^) ' (28) 

i 

^m(^) = E^(^*)>^5g(xj) , (29) 

i.e. the local (L), the cube-smeared (S) and the double cube-smeared (D) axial currents 
|27|. Correlators involving smeared sources were computed in the Coulomb gauge. The 
last two operators were proven to be effective for isolating the lightest 5-meson state. 
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At large time distances the correlation functions above behave as 

CRR\t) ^ Z"^ Z"^' e-^^\ (30) 

and 

In the above equations 

Z^ = -^{0\A§iO,0)\B) , (32) 

and AE is the binding energy of the B-meson. To extract the matrix elements, one 
takes the ratio 

which at large time distances behaves as 



(33) 



where Za is the renormalization constant for the axial current. 
The physical value of the S-parameter is then given by 

Br{m,)= ^ Zo^Zj'Bgf. (35) 

i=L,R,N,S 

We call this method the ratio method for determining the 5-parameter. One may also 
consider the ratio 

i^r E ^o. , (36) 

i=L,R,N,S l*2j 



which at large time distances behaves as 

We call this method the combined-ratio method. 

Both the ratio and the combined-ratio method should give the same value for Bb- 
This is a check of our numerical results. 
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4 Numerical results. 



As explained in the previous section, the determination of Bb requires the computation 
of the two- and three-point correlation functions ( p^ ) and (p7[). The SW-Clover fermion 
action [^] for the light quarks was used, in the quenched approximation. The tree- 
level improved [13, 14 1 propagators of the heavy quarks were computed in the static 
limit. Our results are based on a set of 600 gauge field configurations, computed on a 
lattice of size 24^ x 40 at /3 = 6.0. The calculations were performed at three values of 
the masses of the light quarks, corresponding to K = 0.1425,0.1432 and 0.1440. This 
allows to extrapolate the results to the chiral limit. All the errors have been computed 
with the jacknife method by decimating 30 configurations at a time. 

The procedure to measure Bb is standard. At fixed ti (^2), we study the behaviour 
of the ratios p^ ) and (^) as a function of t2 (ti), searching for a plateau in t2 (ti). 
BBimh) is defined by the weighted average of the data points in the plateau region, if 
this exits. We will take as our best determination of i?B(m;,), the value evaluated in 
a time interval where the ratios appear to be independent of both ti and t2- Notice 
that, contrary to the UKQCD collaboration, we never compute the matrix elements 
of parity-odd terms of the different operators in eqs. (^) and (p^), because they are 
zero by parity. For this reason the matrix element of and Or are equal and in the 
following only the results for Ol will be given. 

In order to improve the isolation of the lightest meson state at short time-distances, 
the ratios (|33|) and (^6|) using single and double cubic smeared axial currents were 
computed. The sizes of the cubes used in our simulation are = 5, 7 and 9. From 
our previous studies, Lg = 7 and 9 were shown to give a good isolation of the lightest 
state for t/a > 4 — 5 or t/a < 36 — 35 [^]. For this reason, even though a plateau 
seems to set in at much ealier times, see fig. |l], the i?-parameters were extracted from 
the ratios (|3^ ) and (|6|) at fixed t2 = 34-36 and for ti > 4. 

In fig. 0, in order to display the quality of our data, the results for i^o^ , Ros 
Ron Rb = Bb as a function of ti/a, at t2 fixed are shown in a specific case. With 
our large statistics, we are able to observe the plateaux over large time-distances. Thus 
we were able to fit the ratios up to ti/a = 7-9. On the basis of our previous experience, 
this makes us confident that the lightest meson state has been isolated. 

We have fitted the ratios ( |3^ ) and ( |36D to a constant, for several time- intervals. As 
an example, in table |2| the value of i?B(mb), at several values of the light-quark Wilson 
parameter K, is given. Since the results are almost independent of the value of the 
light quark mass, they can be safely extrapolated. The extrapolations to the chiral 
limit (Kc = 0.14543(1)) and to the mass of the strange quark (Kg = 0.14367(6)) are 
also presented in table 

Finally, we discuss the dependence of our results on the smearing size. The value 
of all the matrix elements obtained with different smearing sizes are very close. For 
example, in the same case as that considered in table ^, the values of BBirrih) in 
the chiral limit are 0.78(2), 0.76(2) and 0.74(3) for = 5,7,9, respectively. Within 
the statistical errors, a systematic shift of the results with the smearing size cannot 
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Figure 1: The ratios Rol> Ros! Ron ^'^^ from the combined-ratio method are shown as 
a function of the time ti, at t2/a = 35. They were computed at K = 0.1432, by using , 



eq. ^^H), with Ls = 7. The labels SPT and BPT refer to standard and boosted perturbation 
with the V-coupling respectively. The lines show the time-interval of the fits. 



BB{mb) 


K 


t2 = 36, ti = 4-8 


t2 = 35, ti = 4-9 


t2 = 34, ti = 4-6 


0.1425 


0.796(15) 


0.774(19) 


0.79(3) 


0.1432 


0.792(16) 


0.771(19) 


0.79(3) 


0.1440 


0.787(18) 


0.77(2) 


0.79(3) 


0.14543 


0.78(2) 


0.76(2) 


0.79(4) 


0.14367 


0.79(2) 


0.77(2) 


0.79(3) 



Table 2: Values of BB{rnh) obtained as explained in the text, by using , eq. (^), with 
Ls = 7. In order to obtain the renormalized operators, the boosted V-coupling has been used. 
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Figure 2: Comparison of the contributions of the operators Ol, Or, On and Os to the 
B -parameter Bs^iTrib). The renormalized operator has been computed with the V-coupling, 
using the option Mi. 

be appreciated. Thus we added in quadrature the differences of the central values 
obtained with different smearings as an error. In fig. |2|, in order to show the numerical 
importance of the different terms, the contributions of the operators Ol,r,n,s to the 
5-parameter are shown. As expected, the operator Ol gives the largest contribution; 
the correction given by the three other operators is of the order of 20% and cannot be 
neglected. 

We now give the final results by considering first the ratio method. For the different 
operators the results are 

Bo^ = Boj, = (0.94 ± 0.02 ± 0.04) , Bq^, = (0.98 ± 0.05 ± 0.04) , 

Bos = -(0.60 ± 0.01 ± 0.03) , (38) 

where, as usual, the first error is the statistical one and the second is the systematic 
one, which takes into account of differences coming from different choices of the interval 
of the fit and of the smearing size. 

By using in eq. (^) the numerical values of the renormalization constants from 
table |l|, the values of Bb^ given in table I were obtained. The results from the ratio 
or the combined-ratio methods are indistinguishable. Notice, however, the large differ- 
ences coming from different choices of (equivalent) sets of perturbative renormalization 
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Options 


Ratio method 


Combined-ratio method 


Standard 
0.0796 


Ml 


0.87(4) 


0.87(5) 


M2 


0.81(4) 


0.81(5) 


ttg Boosted uq 
0.1458 


Ml 


0.82(5) 


0.81(5) 




0.67(5) 


0.67(4) 


Boosted V 
0.1800 


Ml 


0.77(6) 


0.76(5) 


M2 


0.55(6) 


0.54(4) 



Table 3: Values of BB^ifni,) for different choices of the lattice coupling constants and options, 
see the text. 



constants. From the results of the table, by using eq. (^, it is straightforward to derive 
the renormalization group invariant Bb^- 

In table 0, our results for Bb are compared to other determinations. In the table 

" (1) " (2) 

-B^j and refers to the renormalization group invariant i?-parameters computed 
from BBjirrib) using the leading or next-to-leading formulae respectively. The leading 
formula is obtained from eq. (^ by dropping the term proportional to J5 and by using 
the one-loop expression for Ug (at one loop as(m^) = 0.2342). As for the results of 
the other groups, irrespectively of the fact that they used leading or next-to-leading 
formulae to compute i?B^(mb), we evolved their value from the scale used in the original 
paper to mi, = 5 GeV using the one loop evolution equations. For this reason some of 
the numbers in the table may differ from those quoted in the corresponding references. 

Our results are in good agreement with those of UKQCD pl|, when a similar recipe 
for the renormalization constants is used 0. The UKQCD number reported in table ^ 
should be compared with the value obtained with the M2 recipe and the uq coupling 
constant, see table ^. With the exception of the result of ref. which was obtained 
using the Wilson action for light quarks and is surprisingly large, the static value is 
smaller than other determinations of the same quantity obtained in the "full" theory, 
i.e. by extrapolating in the heavy quark mass mq results obtained for mg < nih |30|- 
|32|. A large value was also obtained by the QCD sum rule calculation of ref. 



^From the study of the S-parameter as a function of the light-quark mass, we also 

get 

Bbs Bb, 



Bb, 



1.01 ± 0.01 . (39) 



5d 

Notice that this ratio is almost independent of choice of the different possible options 
in the calculation of the renormalization constants, see sec. ^. 

Using the same set of gauge field configurations, the following ratio has also been 



^ We recall that there are some differences because they have ommitted the NLO corrections and used by 
mistake the constants of ref. |13[ instead of ref. These differences are, however, numerically unimportant. 
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o(2) 


^BA^h) 


Authors 


Ref. 


Remarks 


1.45(22) 


1.58(24) 


0.99(15) 


Narison et al. 


M 


1 


QSSR 


1.29(7) 


1.40(7) 


0.88(5) 


ELC 


m 


1 


Extrapolated Wilson j3 = 6.4 


1.30(9) 


1.42(10) 


0.89(7) 


Soni et al. 


pi 




Extrapolated Wilson /3 = 5.7 - 6.3 


1.31(7) 


1.42(7) 


0.895(47) 


JLQCD 


ii 




Extrapolated Wilson /? = 6.1 


1.23(9) 


1.34(10) 


0.840(60) 


JLQCD 


p 




Extrapolated Wilson /3 = 6.3 


1.01(6) 


1.10(6) 


0.69(4) 


UKQCD 


ii 




Static Clover j3 = 6.2 


1.42(6) 


1.54(6) 


0.97(4) 


Christensen et al. 


p 


1 


Static Wilson (3 = 6.0 


1.11(7) 


1.21(8) 


0.76(5) 


APE 


This work 


Static Clover (5 = 6.0 Mi 


0.79(6) 


0.86(6) 


0.54(4) 


APE 


This work 


Static Clover (3 = 6.0 M2 



Table 4: Values of the B -parameter as determined by previous studies, and from the present 
one, are presented for comparison. 5^^^ and B^^^ refers to the renormalization group in- 
variant B-parameters computed from Bb^{itli,) using the leading or next-to-leading formulae 
respectively. The results of this work are those obtained by using the "Boosted V" af. 



measured (see ref. [^] for details) 

fB 



1.17 ± 0.03 (40) 
JBa 



Combining eqs. d^) and (II), one gets 

f2 

^= - 1.38 ± 0.07 . (41) 



By using eq. (^ and the above result, one obtains 



(,45 ±0.13)0^. (42) 



where, for the masses and lifetimes of the i3-mesons, we have used the values given in 
ref. l5i. 



5 Conclusions 

In this work we have computed the i?-parameter of the i?-meson in the static limit, 
on 600 gauge field configurations with an improved action. Given the large set of 
configurations, our statistical errors are rather small. The large statistics allows also a 
reduction of the systematic errors coming from the imperfect isolation of the ground 
state. Errors coming from the extrapolation to the chiral limit are also quite small. 
Moreover, for this quantity, the error coming from the quenched approximation was 
estimated to be negligible [^]. Thus the most important source of systematic error 



13 



in our results is the determination of the renormalization constants, which are known 
so far only in one-loop lattice perturbation theory. In order to improve the accuracy 
of the lattice predictions for the B parameter in the static theory, a non-perturbative 
computation of these constants is clearly necessary. 
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